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Abstract. We investigate a simple model for a galactic halo under the assumption 
that it is dominated by a dark matter component in the form of a Bose-Einstein 
condensate involving an ultra-light scalar particle. In particular wc discuss the 
possibility if the dark matter is in superfluid state then a rotating galactic halo might 
contain quantised vortices which would be low-energy analogues of cosmic strings. 
Using known solutions for the density profiles of such vortices we compute the self- 
gravitational interactions in such halos and place bounds on the parameters describing 
such models, such as the mass of the particles involved. 



PACS numbers: 03.75.Nt, 11.27.+d, 95.35.+d 
1. Introduction 

In the standard model of galaxy formation, the visible component of a galaxy is supposed 
to be embedded in an invisible halo of non-baryonic matter [Q [2] . This dark component 
is further supposed to be cold, meaning that it is usually assumed to consist of very heavy 
particles with very low thermal velocities. However, it has been known for some time 
that Cold Dark Matter (CDM) models have certain problems in reproducing observable 
properties of galaxies, among them being the predicted presence of central density cusps 
and the overabundance of small scale structure [3], HI [5]. In the light of these issues, 
some authors (e.g. [6]) have suggested that the Dark Matter could instead consist of 
ultralight particles possessing a de Broglie wavelength sufficiently large that quantum- 
mechanical effects might manifest themselves on astrophysically interesting scales. Such 
models would naturally predict smoother and less centrally concentrated galaxy haloes 
owing than in the CDM case. 

Advocating a particular version of this idea, Silverman & Mallett [7j suggested 
a symmetry breaking mechanism for the production of such a particle, based upon 
a real-valued scalar field. Although in this case the symmetry breaking mechanism 
provides a nice example of particle production in a universe with a cosmological constant, 
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symmetry breaking with a real scalar field generically produces a catastrophic domain 
wall problem [8], and this example would seem to be no exception j9] so this is probably 
not a viable scenario. However, these papers consider the possibility that the Dark 
Matter component resides in a Bose-Einstein Condensate (BEC). The dynamics and 
possible observational consequences of a Cosmological fluid with such properties has 
been investigated [10], using techniques developed in the field of condensed matter 
physics. The equation describing a BEC is known to condensed matter theorists as the 
Gross-Pitaevskii (GP) equation, but is probably more familiar to cosmologists as the 
nonlinear Schrodinger equation (NLSE). 

In condensed matter theory, the term Bose-Einstein Condensate is usually applied 
to a dilute bosonic gas confined by an external potential, the bosons occupying the 
lowest available quantum state. Typically, in the limit of large particle number, the 
density distribution of the condensate is taken to be described by a macroscopic wave- 
function that is considered to be a quantum field. This field is manipulated by the 
Gross-Pitaevskii equation, or nonlinear Schrodinger equation, rather than working with 
the usual creation and annihilation operators of quantum mechanics. The density 
distribution of the condensate can be represented by a macroscopic wave-function of 
the same form as the ground state wave-function of a single particle. The momentum 
distribution of the condensate is obtained by taking the Fourier transform of this wave- 
function. In an experimental setup, the occurrence of a Bose-Einstein condensate is 
confirmed by a sharp peak in the momentum space distribution of the gas of particles. 

More speculatively, the concept of a BEC can also be applied to such hypothetical 
particles as axions or ghosts. In this context, the axion field, for example, is coherent 
and has relatively small spatial gradients. The gradient energy can be interpreted as 
particle momenta, which will be the same and small for each particle, hence giving a 
sharp peak in the momentum space distribution as in the case of the more familiar BEC 
described above. 

In quantum field theory, a condensate corresponds to a non-zero expectation value 
for some operator in the vacuum and, in the limit of large quantum number, this 
condensate can be considered to be a classical field. This is a good model for the 
condensate of Cooper pairs in a superconductor, or for helium atoms in a superfluid 

EH]. 

The usual, linear Schrodinger equation, coupled to the Poisson equation can be 
used to model many phenomena in Cosmology. As well as modelling a quantum 
mechanical system, as in (6], it has also been used as a classical wave equation to 
model structure formation. It has been shown that using the Condensed Matter concept 
of a Madelung transformation to yield the Euler and Continuity equations from the 
Schrodinger equation, applies as well as to a Cosmological fluid as it does to fluids in 
the laboratory [TJl [TSJ El ESI M, E] • 

Silverman & Mallett [7j also considered the rotation of a galactic-scale dark matter 
halo. Using a phenomenological description taken directly from condensed matter, they 
concluded that a galactic halo should be threaded by a lattice of quantised vortices, as 
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a consequence of the rotation of that galaxy. Indeed from studies of rotating BECs and 
quantum turbulence [TBI [19] , it would seem to be difficult to prevent such vortices from 
forming. The galaxy velocity rotation curve produced by these authors reproduces the 
approximate form of observed rotation curves. 

A similar conclusion was reached in Yu and Morgan [20J. This paper considered 
stationary cylindrical solutions of a complex <f) A scalar field model, coupled to gravity. 
These solutions are Nielson-Olesen vortices, also known as local U(l) Cosmic Strings 
[8]. To describe the motion of these vortices in the galaxy, Yu and Morgan's procedure 
was to calculate the motion of one vortex according to a gradient in the phase induced 
by the surrounding vortices. 

There are many models using the Schrodinger-Poisson, or the relativistic Einstein- 
Klein-Gordon, system to describe slightly different physical processes. A non-exhaustive 
list includes scalar field dark matter [2TI [22] . boson stars [23], Oscillatons |24j ; 
condensate stars [25], repulsive dark matter [26] and fluid dark matter [271 [28] . as well 
as the fuzzy dark matter and classical fluid approaches that we have already mentioned, 
and the more established theories such as the Abelian-Higgs model in field theory, and 
the Landau-Ginzberg model in condensed matter. We will not attempt a thorough 
review of each model here, except to say that it is sometimes difficult to explicitly 
distinguish between them. 

The effects of the interaction of gravity with a coherent state of matter, such as a 
BEC, have certainly been considered [291 [30], an d prompted the question of whether it is 
actually possible for DM to be in a coherent quantum state, if the only interaction with 
visible matter is gravitational. Penrose has also used the Schrodinger-Poisson system 
during his 'Quantum State reduction' research program [3~T] . 

In this paper we seek to determine some of the properties of a quantised vortex 
residing in a galactic-scale Bose-Einstein Condensate dark matter. In particular, we 
will place bounds on the parameters that are used to describe such a vortex. For the 
purposes of this paper we presume that the DM does indeed consist of a BEC, formed 
at an earlier stage of Cosmological history and described by the coupled nonlinear 
Schrodinger-Poisson system, and that vortices are present in this cosmological fluid. 

In Section [2] we introduce the basic formalism for describing a BEC using the 
Gross-Pitaevskii (nonlinear Schrodinger) equation, and vortices within it. In Section [3] 
we discuss coupling the NLSE to the Poisson equation. In Sections @] and [5] we look 
at some of the properties of a vortex as a result of gravitational coupling. We present 
some results in Section [6] and a discussion in Section [7J An appendix contains some of 
the approximations we have used in our work, and is referenced in the main body of the 
paper. 

2. Setup 

For our discussion, we use some of the conventions and proceedures set out by Berloff & 
Roberts [32], and Pethick & Smith [TT]. The nonlinear Schrodinger equation is written 
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in the form 

h 2 r 

^ = __V 2 ^ + ^ J \^(x',t)\ 2 V{\x-x'\)dx', (1) 

where m is the mass of a particle in the BEC, and V^(|x — x'\) is the interaction potential 
between bosons. The potential is simplified for a weakly interacting Bose system by 
replacing V(\x — x'\) with a 5-function repulsive potential of strength Vo, giving 

Ti 2 

im t = v 2 m + v m 2 ^. (2) 

2m 

Defining a state that is independent of time to be the 'laboratory frame', \& = 
exp(iE v /h), it is then possible to consider deviations from that state by considering 
the evolution of ip, where x/i = ^ exp(iE v t/h). Here, E v is the chemical potential of 
a boson, in the sense that it is the increase in ground state energy when one boson is 
added to the system. The nonlinear Schrodinger equation used for subsequent analysis 
is then 

iht/H = -|^V 2 ^ + Wl¥ - E v if>. (3) 
2m 

Multiplying equation ([3]) by <ft* and subtracting the complex conjugate of the 
resulting equation we obtain 

d \^\ 2 _ V7 \ 71 , (4) 



dt 



2mi 



We notice that this is of the form of a continuity equation. 



d\^\ 



2 



+ V.(H 2 v)- (5) 

We identify \ip\ 2 as the number density n, and the related momentum density is given 
by 

j = ^(V*v^- W), (6) 

which is equivalent to 

j = mnv. (7) 
This defines for us the mass density, as p = mn = m\^\ 2 , and the velocity 

V 2mi |^| 2 ' 1 ' 

As suggested in the introduction, we can make a 'Madelung transformation' 

ijj = aexp (i<f>J) , (9) 

and, from equation (jSj), we obtain an expression for the velocity of the condensate 

v = -V<A„. (10) 

m 

Here, 0^ is the velocity potential. Substituting the Madelung transformation, and taking 
real and imaginary parts yields the fluid equations: the continuity equation 

«»\VVW=0; (11) 



dt m 
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and the (integrated) Euler equation: 



h 2 V 2 a 1 2 2 . . 

%— = - -mv - V a + E v . (12) 

at 2m a 2 



Often, the identification 



(13) 

V.(a 2 V&/) =0, (14) 



m 

is used, to maintain contact with the more familiar form of the fluid equations 



dt 

%' h 2 V 2 a {V<f>Jf V 2 , E. 



or + (15) 
dt 2m 2 a 2mm 

Here the quantum nature of the fluid is evident only in the first term on the right 

hand side of the second equation, which is often known as the quantum pressure term, 

although dimensionally it is a chemical potential. This term is relevant only on small 

scales, where quantum effects become important, such as in a vortex core, or where the 

condensate meets a boundary. This identification rather hides the quantum nature of 

the fluid with respect to the fluid velocity, which will become particularly relevant when 

we start talking about vortices in the next section. 

By assuming that the condensate reaches a stationary equilibrium state at a 

distance far from any disturbance, equation (jHJ) gives us the relation 

*. = (!)*■ (16) 

When the condensate wave-function reaches a boundary, such as the wall of a container, 
or the core of a vortex is being considered, we can define a distance over which the 
wave-function changes from zero to its bulk value, or where quantum effects become 
important [321 ITT] . 

a = — r (17) 

{2mE v )2 

This is known as the coherence length, or healing length, as it is the distance over which 
the wave-function requires 'healing'. 

2.1. Vortices 

We have already seen that the velocity of the condensate is given by 

v = -V0 W . (18) 
m 

One would expect then, that the condensate would be irrotational, as 

V x (V/) = (19) 

for any scalar, /. This restricts the motion of the condensate much more than a classical 
fluid. The circulation around any contour then, should also be zero. By Stokes' theorm 

T = j>w.d\ = Jjy x v.dA = (20) 
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This condition, known as the Landau state, was first derived in an analysis of superfluid 
Hell [36], and suggests that rotation of such a condensate should be impossible. 
Experiments by Osbourne [37] indicated that the condensate did indeed experience 
rotation. Feynman [38], building on the independent work of Onsager [39], suggested 
that rotation and hence non-zero circulation could be explained by assuming that the 
condensate is threaded by a lattice of parallel vortex lines. It is possible to have 
circulation surrounding a region from which the condensate is excluded, and in this 
case, this would be the vortex core. To see this, we note that the condensate wave- 
function must be single valued, and so around any closed contour, the change in the 
phase of the wave-function A<ft must be a multiple of 2tt. 

A0„ = jv^.dl = 2nl (21) 

where / is an integer. We immediately see that the circulation is quantised in units of 
h/m. 

v.dl = — 2nl = I— (22) 
m m 

To obtain vortex solutions, we work in cylindrical coordinates (r,x,z), and look for 
a static solution of the nonlinear Schrodinger equation, equation To satisfy the 
requirement of single- valuedness, the condensate wave- function must vary as exp(mx), 
with n integer. We make the vortex ansatz 

%l) = R(r) exp(mx). (23) 

It is interesting to note the similarity between this procedure, and that used in obtaining 
Nielson-Olesen vortices, or Cosmic Strings, in the Abelian-Higgs model [8]. This was 
mentioned in Section^ and will be useful shortly for obtaining equation (1271) . as shown in 



Appendix A.l We can obtain an expression for the velocity of a vortex by substituting 



the vortex ansatz (1231) into equation (jSJ) 

fin 1 . . 

v. = X, 24 

r m 

and we note again the discrete nature of the allowed values of velocity. From now on 
we will consider only n = 1 vortices. Vortices with n > 1 are generally expected to be 
unstable, from energy considerations (see for example Chapter 9.2.2 of [II]), and will 
break up into several n — 1 vortices to make up a vortex lattice, as described above. 
We can note further that Cosmic Strings with winding numbers n > 1 are also unstable 
to perturbations [8]. Such defects break down to several n = 1 configurations in both 
a Condensed Matter environment, and a High Energy Field Theoretic one. Feynman 
initially introduced quantised vortices as a purely theoretical tool with which to explain 
the rotation of the condensate, but the experimental verification of the quantisation of 
rotational velocities (e.g. by [40J) demonstrated that these vortices were indeed real. 
The density profile of a vortex (p(r) = m\R(r)\ 2 ) is defined by the vortex equation, 
which results from substituting the vortex ansatz into equation ([3]) 

ft' WR(r) + IdflW 1 i + vs _ = o 

E v 



2mE v 



dr 2 r dr r 5 
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From equation (1161) we see that the density far from the vortex is given by 



Pc 



mRoo = m—. 



(26) 



Analytic solutions of this equation are not known so it must be solved numerically. For 
our anaylses we will use the approximation 



E\ 1/2 

R ( r ) - ( y) i 1 ~ exp(-r/a )] 



(27) 



as discussed in|Appendix A.l 



3. Self-gravity of a BEC Vortex 

In considering Bose-Einstein condensates on scales relevant to structure formation in 
the universe, we must necessarily include gravitational effects. BECs are typically 
sufficiently dilute that the mass densities are not very large, and so a Newtonian 
approximation is sufficient. Gravitational effects can be added to the BEC by including 
a term in the nonlinear Schrodinger equation that couples to the Poisson equation. We 
then have a pair of equations modelling a gravitationally coupled fluid. 

h 2 

ihipt = + V \^\ 2 ip - E v ijj + m<p G ^ (28) 

2m 

V 2 G = 4nGp = 47rGm|^| 2 . (29) 



3.1. Vortices in Gravitationally Coupled BECs 

To obtain vortex solutions, we again work in cylindrical coordinates (r, z), and 
substitute the vortex ansatz if) = R{r)exp{ix) into equations (12"B1) and (12"9~|) . The system 
of equations describing a gravitationally coupled BEC fluid become 

+ ^-R{rf-R(r) + m<p G (r) = 0(30) 
E v 

m) = ^ + ^ = «^W (3D 

Ideally, we would like to find a solution describing the function R(r) in this system, 
so we can compare the density profile of a quantum vortex, to that of one that is 
gravitationally coupled. However, finding a full simultaneous solution to these coupled 
equations is difficult. Firstly, because the nonlinear Schrodinger equation itself is not 
soluble analytically. Secondly, because the vortex density tends to a constant, so the 
Newtonian potential tends to diverge, and thirdly because these equations do not define 
the vortex velocity, which would be providing the centripetal force to withstand the 
gravitational collapse. In other words, all the variables required to provide a fully 
simultaneous static solution are not defined within these two equations. 



2mE v 



d 2 R{r) 
dr 2 



1 dR(r) 

f w 

r dr r 



2 R(r) 
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4. Vortex Stability in Gravitationally Coupled BECs 

Rather than solving the coupled equations (128j) and ( 1291) directly, we can make 
some arguments regarding the stability of a gravitationally coupled BEC vortex, and 
consequently give some bounds on the parameters that describe it. Our analysis is 
based upon consideration of the circular velocity of a BEC vortex, v^r), and the radial 
velocity induced from gravitational attraction, vq(t). t>^(r) is the velocity that the 
vortex density distribution is moving at, for a particular r, while wg(^) would be the 
velocity experienced by a test particle orbiting that density distribution, at a distance r. 
To sustain a vortex, v^ir) must at least be greater than vc(r), otherwise the quantum- 
mechanical forces at work in the vortex are not sufficiently strong to hold itself up 
against gravitational collapse. That is, the vortex is spinning too slowly to provide 
enough centripetal force to balance the gravitational force. For stability, we therefore 
have the bound, 

v u (r) > v G (r) (32) 



4-1- Gravitational Field of a Cylindrically Symmetric System 

To obtain vq(t), we turn to Gauss's Law to determine the gravitational field of a 
cylindrically symmetric mass distribution, and hence obtain the radial gravitational 
velocity of a test particle moving in the field of that system. Gauss's law is 

• dA = -4vrGM cncl (33) 

The density, p(r) = m\R(r)\ 2 , is already determined in terms of the cylindrical r co- 
ordinate, as it is a solution of the vortex equation. The mass enclosed is the density 
pervading a cylinder of radius r and length L. 



M encl = L / 2nrp(r)dr (34) 
Jo 

The left-hand side of Gauss's law, in cylindrical co-ordinates, is 

grdxdz, (35) 



where the integral over the z co-ordinate is again L, the length of the vortex. Gauss's 
law, then, gives us 



gr2nL = -4vrG2vrL / p(r)rdr (36) 
J o 

giving 

g = _^Gm f r lR{r)l 2 rdr _ (37) 
r Jo 

The sign is negative as we have chosen an outward-pointing surface normal in our 
formulation of Gauss's Law, equation (|33|) . which indicates that the gravitational flux 
will always be towards the origin. This leads to the slightly counter-intuitive conclusion 
that a hole (the vortex) in a constant mass density background would seem to produce 
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a gravitational force towards it, but this is really a manifestation of the (extremely) 
thick shell condition. Viewed another way, this static configuration will want to act to 
collapse in, and close the hole. It is this force that is 'unopposed' in equations (1301) and 
(l3Tj) . This need not concern us further, as it is the magnitude that is required for our 
argument. The magnitude of the induced centripetal force is 

9 = ^ (38) 
r 

and the gravitational circular velocity profile vq is given by 

v G {rf = 4nG f p(r)rdr = AixGm f \R{r)\ 2 rdr. (39) 
Jo Jo 

5. Bounds on Parameters 

We now have expressions for vc(r) and v u (r), equations (124)) and (139|) . to go in the 
bound given by equation (J32l) . In Figured] we plot, as an example, v w (r) and vc(r) and 
the density profile for comparison. For this example, we have used values of m = 3.56 
x 10" 59 kg (2 x 1(T 23 eV), E w = 2.5 x 10" 49 J (1.56 x l(r 30 eV) and V = 4.45 x 10" 84 



Jm 3 (3.7 x 10~ 45 eV~ 2 ) as explained in Appendix Appendix A. 2 

The bound on stability, v u (r) > Vc(r), will always be violated at some point, as 
outside the vortex core v u (r) ~ 1/r and fc( r ) ~ r - We must specify what might be an 
acceptable value of r for v u (r) and vs(r) to cross. For a vortex to exist, the density profile 
should be fully established. We take this to mean that the density has essentially reached 



its background level. From the scaled density profile in discussed in Appendix A.l and 
plotted in Figure (1A1|) . we see that the density reaches its background level at a value 
of about ten times the healing length. Using equation ( 1271) in (1391) . equations ( 1391) and 
( 124"|) in fl32|) . and substituting for E v from equation ffTTj) we obtain 

\Z2~7T ( Gfl 2 f n —21 o --2L -2r 2 Tl 



2 \ V a 



2 



2 !_ n 2L _2r 2 —2£- 

2r + 8ra e a o + 8ao e a o — Ira^e a o — ao e a o 



< r--( 4 °) 



mr 



We will fix the healing length ao, and plot Vq against m (fixing ao and m fixes E v , from 
equation (|T7|) ) to give an allowed range of parameter values. We will do this for various 
values of ao, and for various values of r, which we will take to be an integer number of 
healing lengths, r = na^, with the minimum n = 10 as outlined above. Equation (1401) 
then becomes 

Vq > ^Gm 2 n 2 (2n 2 a 2 + Sna 2 e~ n + 8ao 2 e~ n - 2na 2 e~ 2n - a 2 e" 2n ) .(41) 
5.1. Other Bounds 

We can obtain some other bounds to cut off other bits of parameter space. The 
asymptotic vortex density is given by 

*- = -(t)- < 42 > 
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v G (r), v B (r) 
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4.5 
R(r) plotted 
schematically 

4 




x 10 



Figure 1. Velocity Profiles for vg (green) and (blue). Density profile plotted 
schematically for comparison (red). 



If the vortex exists as a component of a galaxy, then there is a minimum and maximum 
density that the vortex can have, given by the maximum and minimum known values 
of mass density within a galaxy: 



Pmin — Poo — Pr 



(43) 

The value of E v in equation ()42|) is fixed (as we are fixing the healing length), and so 
the bound on the density becomes a bound on Vq. 



<v.< 



(44) 



Equation fj4"TT) gives a lower bound on V , so to obtain an upper bound, we use the 
second half of the above relation. 

h 2 



2a 2 p r , 



(45) 
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Another bound is provided because the vortex velocity should never exceed the speed 
of light, 

v u = — <c. (46) 
mr 



It can be seen from equation (1241) that the vortex velocity increases with decreasing 
radius. This relation breaks down within the vortex core, ao, where the vortex velocity 
diverges. Finding an appropriate description is a topic of some interest in Condensed 
Matter theory [33]. We evaluate the maximum vortex velocity at a distance of 5ao from 
the origin, i.e. in a regime where we are sure the relation holds. This gives a bound on 
the mass. 

ft / „N 

m > . (47) 

5cao 

5.2. Values 

To see how the restriction on m and Vq varies, we can think of a range of healing lengths 
that cover all possible scales in a galaxy. 

1 x 10 10 m (3.2 x l(T 10 kpc, ~ 7 x 10~ 2 AU) < a (48) 
a < 1 x 10 22 m (324kpc) (49) 

This range of scales takes us from sub solar system, to that of the largest known galaxies 
(e.g. IC 1101 in the Abell 2029 cluster [M]). At fixed ao we will also cover a large range 
of n; the number of healing lengths where the velocity profiles cross. For the bound 
given in equation (145]) . we take the minimum density found within a galaxy to be the 
cosmological density. This minimum must necessarily be close to the critical density of 
the universe. 

3H 2 

Pmin = Pc = (50) 
07TG 

With H = 70 km s" 1 Mpc" 1 , this gives a value of p m i n = 9.2 x 10 kg m . 



6. Results 

In Figure (j2J), we show a region of the Vq — m parameter space for the healing length 
ao = 1 x 10 16 m (~ 1 pc). The lower bound on Vq is given when v w and vq cross at a 
value of ten times the healing length, n = 10. A vortex could be considered more stable 
if v u and vq cross at a greater number of n, moving us up into the allowed triangular 
region. However, this can soon reach the minimum density bound on Vq. A value of 
n = 10 6 is also plotted, and it is clear that this is outside the bounded region. The lines 
bounding the region of allowed parameter values are given by equations (|4ip . 047p and 

(HSD. 

Figure ([3]) shows allowed regions for various healing lengths, all at a value of n = 10. 
We see that as we move to smaller values of ao, the allowed bounds on m and Vq both 
move up, as expected from equations (l45l) and ( )47|) . More physically, as the mass of the 
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(eV" 2 ) (eVa Q 3 ) (Jm 3 ) 

8.2e-22 6.2e-90 1e-60 - n=1e6 

8.2e-27 6.2e-95 1e-65 - 

8.2e-32 6.2e-100 1e-70 - 

8.2e-37 6.2e-105 1e-75 - 



8.2e-42 6.2e-110 1e-80 - 

8.2e-47 6.2e-115 1e-85 - 

8.2e-52 6.2e-120 1e-90 - 

8.2e-57 6.2e-125 1e-95 - 

8.2e-62 6.2e-130 1e-100 - 

8.2e-67 6.2e-135 1e-105 - 

1e-62 1e-60 1e-58 1e-56 1e-54 1e-52 1e-50 1e-48 (kg) 
5.6e-27 5.6e-25 5.6e-23 5.6e-21 5.6e-19 5.6e-17 5.6e-15 5.6e-13 (eV) 

Figure 2. Allowed region in Vq — m parameter space, for a healing length of 
a = 1 x 10 16 m (~ 1 parsec) 

particle is increased, the repulsive potential Vq must increase to balance the stronger 
gravitational force. 

7. Discussion 

In this paper we have used techniques from condensed matter theory in a cosmological 
setting to place bounds on parameters describing a Dark Matter candidate, on the 
assumption that the Dark Matter halo consists of a Bose-Einstein Condensate, in which 
quantised vortices reside. In the case of a laboratory BEC, self-gravitational forces 
are not important and even in that case analytical progress is limited. Using a simple 
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(eV~ 2 ) (Jm 3 ) 



8.2e-27 1e-65 



8.2e-32 1e-70 



8.2e-37 1e-75 



8.2e-42 1e-80 



8.2e-47 1e-85 



8.2e-52 1e-90 



8.2e-57 1e-95 



8.2e-621e-100 




1e22 m 



I I 



1 e-65 1 e-63 1 e-61 1 e-59 1 e-57 1 e-55 1 e-53 1 e-51 1 e-49 1 e-47 1 e-45 1 e-43 1 e-41 1 e-39 1 e-37 (kg) 
5.6e-30 5.6e-26 5.6e-22 5.6e-18 5.6e-14 5.6e-10 5.6e-6 5.6e-2 (eV) 



Figure 3. Allowed regions in Vq — m parameter space, with n = 10. Healing lengths 
as labelled. 



physical argument, however, we have shown how rough limits on the consistency of such 
a model can be imposed. Considering a Dark Matter particle of a particular mass, and 
a vortex of a certain radius, places constraints on the values that the chemical potential, 
and interaction potential can take. There remain sizeable regions of parameter space in 
which the model appears to be viable. 

In future work, it would be interesting to investigate further whether a Dark 
Matter candidate could reside in a coherent quantum state, if the only interaction 
was gravitational. A less ambitious undertaking would be to see if the Madelung 
transformation provides a solution to the problem of defining all the relevant variables, 
as suggested in Section 13.11 This would give a set of fluid equations that includes the 
velocity giving rise to the stabilising centripetal force. One problem to be anticipated in 
such a solution, would be that the velocity in the vortex core would still be ill-defined, 
as alluded to in Section El The system would therefore have to be solved by a more 
complete numerical method than we have been able to implement so far. 
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Appendix A. Approximations 

Appendix A.l. Approximations to the Density Profile 

The numerical solution to the NLSE can be cumbersome to work with, so we provide 
some discussion of some approximations that can be used. It is possible to scale the 
the variables r and R(r) in equation (|25l) to obtain a scale-free equation. Scaling r by 
the healing length, r' = r/a , and R(r) by the steady state value, R'{r') = i?(r)/i?oo we 
obtain 

Our first idea for an approximation comes from the field of cosmic strings. The method 
of approximation is detailed in [5] . Looking at the profile of the Higgs field in a Nielson- 
Olesen vortex we see that it can be written, in a similarly scaled way, as 

^ + - j^/MWO - 1) 2 - \fV)UV? - 1) = o ( a.2) 

Here a is a gauge term arising from the coupling to Electromagnetism, and A is 
determined by the potential term of the theory. It is possible to linearise equation 
(1A.2|) to obtain a modified Bessel function as the first order approximation to f'( r ') - the 



zeroth order being 1. This happens in the string case, because the gauge contributions 
serve to cancel one of the terms, leaving the modified Bessel's equation. The linearised 
version of equation ( 1A.1I) does not quite reduce to a modified Bessel's equation, but 



taking our lead from the cosmic string example, we write 

R'{r') ~ 1 -exp(-r'). (A.3) 

Another approximation, which might seem to be more accurate, was developed 
by Berloff [35] in a condensed matter context. The Pade approximation has the same 
asymptotics at r = and r = oo as the function one is trying to approximate. The 
Pade approximation in this case gives 



r 



/2 



(0.3437 + 0.0286r' 2 ) 



R'(r')~\ v I V (A.4) 

v ' V 1 + 0.3333r /2 + 0.0286r' 4 v ; 

This solution is plotted in Figure IA~T1 along with the numeric solution given by equation 

(lA.lj) . and the previous approximation, equation flA.3j) . The Pade approximation is 



indeed much more accurate in the small and large r regions. However, the Pade 
approximation has the tendency to overestimate the density in the central region, 
producing a density function whose derivative is negative in this region. As discussed in 
the main body of this paper, the gravitational potential is proportional to the density, 
and so the gravitational force will be proportional to the derivative of the density 
function. If we chose to use the Pade approximation for our density profile, we could be 
potentially misled by its behaviour in the central region. 
We will use the approximation 



i 



R(r)= ('^) 2 [l-exp[-r/a ]). (A.5) 
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1.2 

R'(r') 




01 23456789 10 

r' 



Figure Al. Numeric solution to Equation (jA.lj) (blue), the Pade approximation 
equation (|A.4j) (red), and the scaled approximation used in this analysis, equation 
(|A.3D (green). 

Appendix A. 2. Approximations for Parameters Defining the BEC 

To enable us to obtain actual values for the velocity and density profiles that we are 
considering, we must provide values for the parameters m, Vq, and E v . The properties 
of Dark Matter particles are, by their very nature, unknown, so we must make some 
approximations. We use the analysis in [7j to provide us with some data values. The 
mass of the Bose Einstein Condensate Dark Matter particle in that paper is 3.56 x 
10~ 59 kg (2 x 10~ 23 eV). Their analysis is based on the mass and angular rotation 
of the Andromeda galaxy. The mean density is given as 2 x 10 _24 kg m -3 , and they 
estimate that the vortex line density in the galaxy would be about 1 vortex per 208 kpc 2 . 
This gives a vortex radius of r w ~ 2.5 x 10 20 m. We again turn to vortex lattices in 
condensed matter systems to provide us with some further estimates of vortex properties 
in a BEC. 

Taking the distance between two vortices to be twice the vortex radius, we note 
from experimental observations of vortex lattices in a BEC that the vortex density 
reaches the normal density at about half the vortex radius; see, for example, Figure 9.3 
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in [TT], taken from [4"T| . From Figure (1A1I) . we also see that the vortex density reaches 
the normal condensate density at around five healing lengths. This gives us an estimate 
of = 5ao- We then use ~ 2.5 x 10 20 m, ao = %/ {2mE v )^ , and poo = mE v /Vo to 
give estimates for and V . With these approximations we find values of E v = 2.5 x 
1(T 49 J (1.56 x l(T 30 eV) and V = 4.45x 10" 84 J m 3 (3.7 x 10" 45 eV" 2 ). 
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